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We provide integrated cross sections for quasielastic charged-current neutrino-nucleus scattering. 
Results evaluated using the phenomenological scaling function extracted from the analysis of ex- 
perimental (e, e') data are compared with those obtained within the framework of the relativistic 
impulse approximation. We show that very reasonable agreement is reached when a description 
of final-state interactions based on the relativistic mean field is included. This is consistent with 
previous studies of differential cross sections which are in accord with the universality property of 
the superscaling function. 
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The development of present and future experimental 
studies of neutrino oscillations at intermediate to high en- 
ergies benefits from high-quality predictions for neutrino- 
nucleus cross sections. The kinematics involved in these 
processes typically lie in a domain where a fully relativis- 
tic formalism is required; not only should the reaction 
mechanism incorporate relativity, but also the nuclear 
dynamics must be described in a relativistic framework. 
Moreover, any reliable model of neutrino-nucleus scat- 
tering first needs to be tested against electron scatter- 
ing for similar kinematics. Fully relativistic results for 
the latter have been presented in the past p], [H, Q us- 
ing the impulse approximation, namely, where only one 
nucleon in the nucleus interacts with the virtual photon 
exchanged in the process. Such analyses have been shown 
to provide accurate descriptions of quasielastic (QE) ex- 
clusive (e, e'p) processes when a proper description of the 
final-state interactions (FSI) between the ejected nucleon 
and the residual nucleus is incorporated. This is accom- 
plished by using complex relativistic optical potentials 
which have been fitted to elastic nucleon-nucleus scatter- 
ing data For inclusive processes of the type (e, e') and 
{y, (a), the contribution of all channels must be retained; 
hence the use of complex potentials should be avoided be- 
cause of the loss of flux implied by the imaginary term. 
Different approaches have been considered in the liter- 
ature, e.g., the use of purely real potentials within the 
impulse approximation 0, (g,J3 and analyses based on 
the Green function method 

[EH G3- Although both 
treatments lead to similar results, a comparison with ex- 
perimental QE cross sections is not yet conclusive be- 



cause of the effects introduced by ingredients beyond the 
impulse approximation such as long- and short-range cor- 
relations, meson-exchange currents (MEC) and, at high 
energies, the excitation of the A resonance. 

These difficulties can be partially overcome by making 
use of the scaling behavior of the electron-nucleus cross 
sections through a SuperScaling Analysis (SuSA) in the 
region of both the QE and A peaks. Previous investiga- 
tions of inclusive (e, e') world data have clearly demon- 
strated the validity of scaling and sup erscaling properties 
in these kinematical domains [ll|, [l2|, HH ■ To summarize: 
by dividing the experimental (e, e') differential cross sec- 
tions by an appropriate single-nucleon factor one gets 
the superscaling function, which embodies the basic in- 
formation about the nuclear dynamics. At sufficiently 
high energies this function depends upon the transferred 
momentum (q) and energy (to) only through a particular 
combination, the scaling variable Tp(q,uj) [H, (first- 
kind scaling) and is independent of the particular nucleus 
selected (second- kind scaling). 

Importantly, the superscaling function extracted from 
data presents an asymmetric shape with a pronounced 
tail extending into the region of high transferred ener- 
gies, corresponding to positive values of the scaling vari- 
able tp. While this asymmetry is largely absent in most 
non-relativistic models based on the impulse approxima- 
tion, the systematic study presented in [f| @] has shown 
that the correct amount of asymmetry is provided in the 
relativistic impulse approximation (RIA) when FSI are 
described with a relativistic mean field (RMF) potential. 
Recently, a similar asymmetric scaling function has also 



2 



been obtained within a semi-relativistic (SR) approach 
including FSI through a Dirac-equation-based model [3] • 

The superscaling function /(V>) has been investigated 
in detail using various models based on the impulse ap- 
proximation. The results are consistent with the exis- 
tence of a universal scaling function, which means that 
/(VO is basically the same for different types of reactions, 
such as (e, e') and (y, /x), provided that the same kine- 
matical regime is considered. As a consequence, instead 
of using specific models for the nuclear structure and 
reaction mechanism, one can use the phenomenological 
SuSA superscaling function extracted from (e, e') data to 
make reliable predictions for neutrino-nucleus cross sec- 
tions. This strategy was first pursued in [l5| for charge- 
changing (CC) (v, ji) processes, and has recently been ap- 
plied [16| to neutral-current (NC) neutrino-nucleus scat- 
tering reactions. In both cases, predictions were given for 
differential cross sections and a comparison with various 
relativistic FSI predictions was also provided for the two 
i-channel processes, (e, e') and (v,fi). 

The differential cross sections and scaling functions 
were studied at depth in [f| 0, HE E3] , while in this work 
we focus on integrated cross sections. We present and 
discuss results for CC neutrino cross sections, integrated 
over the muon energy as a function of the muon scat- 
tering angle, as well as for cross sections integrated with 
respect to the scattering angle as a function of the out- 
going muon energy. Finally, the behavior of the total 
cross section (integrated over both the scattering angle 
and muon energy) as a function of the incident neutrino 
energy is also investigated. We compare the results eval- 
uated within the RIA framework with those obtained by 
using the phenomenological scaling function. We also in- 
clude for reference the results corresponding to the rela- 
tivistic Fermi gas (RFG) model. 

Studies of integrated cross sections are relevant for the 
analysis of present and future neutrino oscillation exper- 
iments. The significant differences observed in the differ- 
ential cross sections when comparing various model pre- 
dictions [1, 0] with the phenomenological scaling function 
can also have consequences for the integrated cross sec- 
tions. In particular, it is of interest to determine how 
the large asymmetry in the scaling function implied by 
the SuSA approach reflects on the integrated cross sec- 
tion. A comparison with the different theoretical descrip- 
tions, specifically the RMF approach which reproduces 
the asymmetry of data, may give us important clues on 
the validity of the diverse models. 

In performing the analysis of integrated cross sections, 
an important issue refers to some of the kinematical re- 
gions accessed experimentally. In particular, small values 
of muon scattering angle between the incoming neutrino 
and the outgoing muon, 6^, imply small values of the mo- 
mentum and energy transferred in the process. In this 
situation, a theoretical description of the reaction mech- 
anism based on the impulse approximation is question- 
able. Fortunately, the contribution at these small angles 
to the global integrated cross section is negligible because 



of the phase-space factor sin 9^ entering the integral. Our 
aim here is to estimate the effects introduced in the inte- 
grated cross sections by different theoretical models and 
to compare these with those obtained by using the SuSA 
approach, and thus we integrate over the full range of (9 M 
allowed by the kinematics. 
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FIG. 1: CC cross section for 12 C integrated over the muon 
scattering angle as a function of the kinetic energy of the 
emitted muon. Results for RMF (solid), RFG (dashed), SuSA 
(dotted), and SRWS (dot-dashed) are compared for three dif- 
ferent choices of kinematics. 



In Fig. Q] we present results for the CC 12 C(f M ,^~) 
cross section integrated over the muon scattering angle 
as a function of the muon kinetic energy. Predictions 
corresponding to the relativistic impulse approximation 
with FSI described with the RMF potential (solid line) 
are compared with the cross sections evaluated by making 
use of the phenomenological superscaling function f(ip) 
(dotted line), labelled SuSA. We also show results ob- 
tained in the RFG model (dashed line) and in the SR 
approach described in [l~7| using a Woods Saxon potential 
(dot-dashed line) denoted as SRWS. The parameters of 
the WS potential have been taken from [I?) , but with the 
depth of the central part of the proton potential adjusted 
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to reproduce the experimental mass difference between 
12 C and 12 N. In the RMF and SRWS approaches the 
experimental energies for the bound neutrons are used, 
whereas in the SuSA and RFG approaches we have taken 
into account the mass difference between the initial and 
final nuclear systems, in accord with with Eq. (6) in [H[. 
This cuts out the higher muon kinetic energies and pro- 
duces a depletion of the cross sections. 

Three different values of the incident neutrino energy 
have been considered: £„ = 0.5, 1 and 1.5 GeV. The rel- 
ativistic plane-wave limit for the final nucleons, i.e., no 
FSI considered (not shown in the figure), leads to cross 
sections which are close to the RFG curves. Results in 
Fig. [1] show that the inclusion of FSI in SuSA and RMF 
gives rise to a significant shift of strength: the cross sec- 
tion increases the RFG result at lower muon energy val- 
ues, whereas it quenches it at higher energies. This re- 
duction, similar for both RMF and SuSA approaches, is 
about 20-25% in the region close to the maximum. In 
contrast, the SR model gives results that are similar to 
those of the RFG. Note also that the use of real poten- 
tials for describing the final nucleon states leads to the 
resonant structure observed for high T M (that is, small 
energy transfer ui). 

In Fig. [5] we present the cross section integrated over 
the muon kinetic energy as a function of the muon scat- 
tering angle. Again we compare results for RMF (solid 
line), RFG (dashed), SuSA (dotted) and SRWS (dot- 
dashed). The inclusion of FSI in RMF and SuSA gives 
rise to a depletion of the cross section, except for very 
small angles. 

As we show in the following, this implies a smaller 
global integrated cross section, which means that, be- 
sides the significant redistribution of strength produced 
by FSI, a global reduction is also observed when an in- 
tegral over the muon energy and scattering angle is per- 
formed. 

From Figs. [T]and [2] it appears that the SuSA and RMF 
approaches yield very similar results, corresponding to a 
redistribution of strength with respect to the RFG and 
SR models, in agreement with what was found in Q and 
pjj for the double-differential cross sections. Note also 
that nuclear model effects are more important for lower 
neutrino energies (e„ = 0.5 GeV), in accord with previous 
work [EEl], and that the scaling arguments are only valid 
for high enough values of the transfer momentum q, i.e., 
for high enough values of the neutrino energy. Hence, 
some caution must be exercised when using the SuSA 
approach for small s u . 

Finally, nuclear model and FSI effects in the fully in- 
tegrated cross section are presented in Fig. [3l where the 
total cross section a is plotted as a function of the inci- 
dent neutrino energy. Besides the four models considered 
in the previous figures we also display the result corre- 
sponding to RPWIA. All of the descriptions lead to a 
similar behavior for the quasielastic cross section, which 
increases with the neutrino energy up to e v ~ 1-1.2 GeV 
and then saturates to an almost constant value [H, [l|| . 
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FIG. 2: Cross sections for the reaction 12 C(^ M ,// _ ) inte- 
grated over the emitted muon energy as a function of the 
scattering angle 6^ for three values of the incident neutrino 
energy e„. Results obtained within the context of the RIA 
and FSI described through the RMF (solid lines) and SRWS 
(dot-dashed) approaches are compared with the RFG model 
(dashed) and with the cross section obtained by using the 
phenomenological superscaling function extracted from the 
analysis of (e, e') world data, denoted by SuSA (dotted line). 



We observe that the RPWIA and RFG models give 
very similar results up to e v ~ 0.8-1 GeV, while for 
higher energies the RFG yields a lower cross section due 
to the above-mentioned nuclear mass difference, which 
entails a cut in the lower energy transfers. The SRWS 
model yields results which are close to the RFG predic- 
tion. The results denoted SRWS-tot include also the con- 
tribution of the discrete spectrum of 12 N obtained with 
the WS potential. This contribution turns out to be very 
small (below 2%). On the other hand, the RMF predic- 
tion coincides with the SuSA one up to neutrino energies 
of about 1.2 GeV. Moreover, the FSI effects in RMF and 
SuSA remain sizeable even at large £„. In fact, they tend 
to stabilize, being of the order of ~15% at neutrino en- 
ergies e v > 1.2 GeV. 
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FIG. 3: (Color online) Total integrated CC cross section 
a for QE muon neutrino reactions with 12 C as a function 
of the incident neutrino energy. We present results corre- 
sponding to RMF (squares), RFG (solid line), SuSA (dashed 
line), RPWIA (dot-dashed line), SRWS (circles) and SRWS- 
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Before concluding, a comment is in order concerning 
the effect of Pauli blocking (PB) in the RFG model. It 
is well-known that PB, which is obviously accounted for 
in the RFG model, only affects the low momentum and 
energy transfer region. However in the process we are 
considering here this region is kinematically forbidden 
due to the mass difference between the initial ( 12 C) and 
residual ( 12 N) nuclei. As a consequence the effect of PB 
in the integrated cross sections turns out to be negligible. 

Although no direct comparison of our predictions with 
experiment can be performed, since no data on 12 C are 
yet available in this energy range, the existing experimen- 
tal data on 2 H (which are, e.g., summarized in Ref. (20j ) 
are in qualitative agreement with the SuSA predictions. 
However, due to the above mentioned relevance of the 
mass difference bewteen the initial and final nucleus, a 
simple renormalization of the nucleon number is not suf- 



ficient to test our model based on deuteron data. 

In conclusion, we have evaluated the charged-current 
quasielastic neutrino-nucleus cross sections integrated 
over the muon scattering variables (kinetic energy and 
scattering angle) within different relativistic theoretical 
approaches. Our results can be summarized as follows: 
(i) The effect of nuclear interactions is sizable for all val- 
ues of the neutrino energy ranging from 0.5 to 2 GeV and 
amounts to a significant redistribution of the strength in 
the single-differential cross sections and to a lowering by 
about 15-20% of the total cross section with respect to 
the relativistic Fermi gas result, (ii) At high neutrino 
energies the differential cross section da/dO^ is strongly 
peaked at low scattering angles; similarly da/dT^ dis- 
plays a pronounced maximum at high muon kinetic en- 
ergy. As a consequence the total cross section is rather 
sensitive to the range of integration, which should be 
carefully taken into account when comparing with ex- 
perimental data, (iii) Finally, the integrated cross sec- 
tion evaluated with the phenomenological superscaling 
function is very close to the RMF prediction. This result 
complements our previous analyses of differential cross 
sections and scaling functions, and it gives us confidence 
in the adequacy of descriptions of QE {y ^ , /U reactions 
within the RIA context when FSI arc included via strong 
relativistic potentials. 
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